Energy Conversion and Management 79 (2014) 666-673 



ELSEVIER 


Contents lists available at ScienceDirect 

Energy Conversion and Management 

journal homepage: www.elsevier.com/locate/enconman 



Thermodynamic model to study a solar collector for its application 
to Stirling engines 

Amir Abdollahpour , Mohammad H. Ahmadi • , Amir H. Mohammadi • • 

a Faculty of Mechanical Engineering, K.N. Toosi University, Tehran, Iran 

b Department of Renewable Energies, Faculty of New Science and Technologies, University of Tehran, Tehran, Iran 
c Institut de Recherche en Genie Chimique et Petrolier (IRGCP), Paris Cedex, France 
d Thermodynamics Research Unit, School of Chemical Engineering, University of KwaZulu-Natal, Howard College Campus, King George V Avenue, Durban 4041, South Africa 


& 


> CrossMark 


ARTICLE INFO 


ABSTRACT 


Article history: 

Received 22 August 2013 
Accepted 16 December 2013 
Available online 31 January 2014 


Keywords: 

Solar collector 
Optical analysis 
Overall loss coefficient 
Fluid temperature 
Absorber tube 
Heat loss 
Thermal analysis 


Energy production through clean and green sources has been paid attention over the last decades owing 
to high energy consumption and environmental emission. Solar energy is one of the most useful energy 
sources. Due to high investment cost of centralized generation of electricity and considerable loss in the 
network, it is necessary to look forward to decentralized electricity generation technologies. Stirling 
engines have high efficiency and are able to be coupled with solar energy which cannot be applied in 
internal combustion engines. Solar Stirling engines can be commercialized and used to generate decen¬ 
tralized electricity in small to medium levels. One of the most important steps to set up an efficient solar 
Stirling engine is choosing and designing the collector. In this study, a solar parabolic collector with 
3500 W of power for its application to Stirling engines was designed and analyzed (It is the thermal inlet 
power for a Stirling engine). We studied the parabolic collector based on optical and thermal analysis. In 
this case, solar energy is focused by a concentrating mirror and transferred to a pipe containing fluid. 
MATLAB software was used for obtaining the parameters of the collector, with respect to the geographic, 
temporal, and environmental conditions, fluid inlet temperature and some other considerations. After 
obtaining the results of the design, we studied the effects of changing some conditions and parameters 
such as annular space pressure, type of the gas, wind velocity, environment temperature and absorber 
pipe coating. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

As a non-emission and lasting energy, solar energy and its 
related applications has been widely studied such as solar build¬ 
ings [1], solar water heating systems and solar energy conversion 
systems [2-19 . In solar energy systems, thermal energy storage 
(TES) unit is an important component to ensure the system stabil¬ 
ity and steady operation with high efficiency performance; due to 
dependence of its performance on weather conditions and time. 
Latent heat thermal energy storage (LHTES) has more density 
and grater temperature of constant phase change, in comparison 
with TES. Recently, parabolic trough power plant has become as 
a unique and progressed mean for solar electricity production. A 
parabolic trough solar collector (PTC) utilizes the energy of sun 
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radiation and employs its useful thermal energy for the heat trans¬ 
fer fluid (HTF) which circulates in the system. 

By determining the geometry and thermal properties of a solar 
system, the thermal performance and energy obtained by the HTF 
can be estimated under different conditions. PTCs usually operate 
at 400 °C. Synthetic oil is utilized as HTF generally. Primary 
design of the solar power plant, prediction of thermal losses, study 
of collector degradation effects and HTF flow rate control require 
heat transfer analysis [20 . Researchers have been working on 
numerous heat transfer models of PTCs, since 1970s. 

Edenburn [21,22] evaluated the efficiency of a PTC by using an 
analytical heat transfer method and the results were satisfying. 
Ratzel et al. [23] performed both analytical and numerical study 
of the heat losses in an annular receiver. Clark [24] has studied 
the thermal and economical performance of parabolic trough 
receivers. Dudley et al. [25] presented a one dimensional steady 
state thermal model of SEGS LS-2 parabolic collector. This model 
was verified with experimental data of Sandia National Laborato¬ 
ries (SNL) for various receiver annulus settings. 

Thomas [26] has investigated a set of numerical equations 
based on heat transfer and heat loss analysis in a PTC receiver. A 
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Nomenclature 


Af the ratio of loss area to total area of concentrator 

C p specific heat capacity (J kg -1 °C _1 ) 

D diameter (m) 

F r the ratio of heat transfer to maximum heat transfer 

F collector efficiency factor 

h convection heat transfer coefficient (W m -2 °C _1 ) 

hf HTF convection heat transfer coefficient at 7) (W nrr 2 °- 

C" 1 ) 

I b beam radiation incident on the receiver surface 

k thermal conductivity (W m -1 °C _1 ) 

g gravitational constant (9.81 m s~ 2 ) 

m mass flow rate (kg s -1 ) 

Nu Nusselt number 

Pr Prandtl number 

Ra Rayleigh number 

T temperature (°C) 

T p is the mean temperature of absorber tube 

N number of day from 1st January 

h s the hour angle 

La latitude 

U L achieved through heat transfer modeling between 

external surface of absorber tube and environment 


W a the collector aperture (m) 

Greek letters 

v kinematic viscosity (m 2 s _1 ) 

P volumetric thermal expansion coefficient (ideal gas) 

(IC 1 ) 

g Stefan-Boltzmann constant 

(p r rim angle 

(p the angle between collector axis and reflector beam 

p surface reflectance 

g emissivity 

a absorbance 

8 S inclination angle 

Acronym 

SNL Sandia National Laboratories 

TES thermal energy storage 

LHTES latent heat thermal energy storage 
PTC parabolic trough solar collector 

HTF heat transfer fluid 


one dimensional energy balanced model was designed by Forristall 
to evaluate heat transfer of solar receiver [27 . Stuetzle [28] 
presented an unsteady state analysis of solar collector by partial 
differential equations to determine the outlet temperature. Vallad- 
ares and Velasquez 29] investigated a numerical model for a single 
pass solar receiver. They showed that the proposed configuration 
increases the thermal efficiency of the solar collector. Recently, 
three dimensional heat transfer analysis of PTCs was carried out 
by simultaneously use of the Monte Carlo Ray Trace Method 
(MCRT) and CFD analysis [11,30,31]. 

In the present work, initially, we review the solar radiation 
energy and the effective parameters on absorbed solar flux. These 
parameters depend on season, time of day, location and incident 
beam. Then, we analyze based on optical. Absorbed reflective flux 
by absorber pipe is one of the results of the optical analysis. In the 
next step, thermal analysis with respect to the energy equilibrium 
model is carried out. This model includes all the necessary equa¬ 
tions to determine the components in the energy balance such as 
type of the collector, terms of absorber pipe, the optical properties 
and the environmental conditions. The absorbed heat is results of 
this analysis that are gained using MATLAB software. According 
to the analysis and the considered power of collector, we wrote a 
code in MATLAB software and we studied and designed the collec¬ 
tor. This design is considered the geographic, temporal, and envi¬ 
ronmental conditions such as the atmospheric pressure, the 
environmental temperature, and wind velocity. Also, other factors 
like type of the fluid, the annular space pressure (between the ab¬ 
sorber pipe and the glass coating) and some design considerations 
for the components of the collector including concentrator, absor¬ 
ber tube and glass coating on the absorber are studied and 
designed. 


2. Optical analysis 

Kalogirou studied parabolic collectors through optical analysis. 
The angle between incident beam on the edge of concentrating 
surface (where mirror radius r r , is maximum) and center line of 
the collector is cp r which is called rim angle [32]. 


Diameter is found from following equation. 
D = 2r r sin(0 m ) 



6 m which is shown in Fig. 1, is a function of the tracking mechanism 
precision and reflective surface disorders [32]. 

r =- — - (2) 

1 + cos (cp) y ' 

where cp is the angle between the collector axis and the reflector 
beam. 

When (p changes from 0 to cp r consequently, r changes from/to 
r r and D increases from 2/sin(0 m ) to £q. (2) in (p r angle; 


r r 


2 / 

1 + cos (cp r ) 


W a is found from Eq. (4) 



W a = 2r r sin(cp r ) (4) 

where W a is the collector aperture (m). 

Through Eqs. (3) and (4) 



4/ sin(<jo r ) 

1 + cos (cp r ) 



W a = 4/tan (^t) 



In parabolic collectors, a number of reflective beams to the 
concentrating do not have incident to receiver, which is known 
as end effect and the area is illustrated in Fig. 2. 


A e = JW a tan(6) 


W 2 

1 + a 

^ 48/ 2 



These kinds of collectors are often restricted by opaque plates 
which causes shadow on a part of reflector [32,38,39 . 



2 

3 


W a h p tan(6) 
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The total loss area is indicated in the following equation 

[32,38,39]: 

A l =(^W a +JW a ( 1+^)) (9) 

Af is the ratio of loss area to the total area of concentrator. The effi¬ 
ciency increases with a decrease in the ratio. 

Af= x a (10) 


2.2. Absorbed reflective flux by absorber tube 

Sukhatme and Nayak conducted a thermal analysis research on 
collectors [33 . All incident beams are not absorbed by the absor¬ 
ber tube. In Eq. (11) S is the amount of absorbed flux by the absor¬ 
ber tube [32,38,39]. 

s = hpRby(m) b + I b R b (m) b ( w f°_ D ) (11) 


3. Thermal analysis 


The efficient absorbed heat is achieved using the energy conver¬ 
sion equilibrium for the inlet and outlet temperatures 

Qu = rhfC p (Tf j0Ut - Tf j n ) (13) 

where T fmin and T f}0Ut are inlet and outlet fluid temperatures, respec¬ 
tively. C p and rhf are thermal capacity and mass flow rate, respec¬ 
tively. With respect to absorber thermal loss [32,38]; 


Qu = F R (W a 

C _ Wq-Dp 

tcD 0 

rh/Cp 

* tcD 0 LU l 


D 0 )L[s - Ui(Tf j n - T a )/C] 


1 



FnD 0 LU L \\ 
m f C p )) 


(14) 

(15) 

(16) 


F 



1 

_1_i_ Do 

U L ^ Djh f 



F r is the ratio of heat transfer to the maximum heat transfer. The 
steady state energy equilibrium for the absorber tube is shown in 
the following equation: 


Qu = (W a - D 0 )Ls - U l tcD 0 L(T p - T a )\ (18) 

U L is achieved through heat transfer modeling between the external 
surface of the absorber pipe and environment, T p is the mean tem¬ 
perature of absorber pipe and D 0 is the tube outside diameter (m). 
The heat transfer is steady and one dimensional, assuming that 
the absorber and cover are in an infinitely long concentric position. 


3.2. Energy equilibrium model 

The steady state one-dimensional energy balance for a point of 
the absorber is illustrated in Fig. 3a and thermal resistance model 
is shown in Fig. 3b. 


where the first part is the amount of reflected absorbed radiation by 

the pipe and the second one is the absorbed solar radiation. 3Z Convention heat transfer between absorber tube and fluid 

R b is indicated in Eq. (12) 

Through Newton’s law of cooling 

R b - . , . s C0S ] -J- J- (12) Quconv = h,D 2 n(T 2 - Tj) (19) 

sin La sin S s + cos La cos S s cos h s 


3 S is inclination angle and obtained from this equation 

[32,33,39,40]: 

sin d s = 23.45 cos[0.986(N + 284)] (12a) 


hi = Nu D2 (20) 

E>2 

Nusselt number is a function of the type of the flow and in the 
operating conditions; the flow is turbulent [33]. 


N is number of days from 1st January, h s is hour angle. This angle is 
based on rotational motion of the Earth in 24 h. Thus, 1 h is equal to 
15° and La is latitude. 


3.2.2. Transient and turbulent flow 

Nusselt number for transient and turbulent flow can be found 
with the following equations [33] 



f 2 /8(Re D2 - 1000)Pr! 

/PrA 

1+12.7^/2/81 

fPrf-l) 

vv 


/ 2 = (1.821og 10 (Re D2 )-1.64)- 2 


( 21 ) 

( 22 ) 


3.2.2. Laminar flow 

In laminar flow, when Reynolds number is less than 2300, Nus¬ 
selt number is 4.36 and constant [34 . 
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With Glass Envelope Without Glass Envelope 



radiation 


radiation 


(l) -AAA/V- ( 2 ) -A/VW 0 ) (4) -A/WV (5) 

convection conduction v v v . conduction 

convection 


xA/VVV* 

wvw 

convection 


(7) 

( 6 ) 


(1) heat transfer fluid 

(2) absorber inner surface 

(3) absorber outer surface 

(4) glass envelope inner surface 


(5) glass envelope outer surface 

(6) surrounding air 

(7) sky 


Fig. 3. (a and b) Steady state one-dimensional energy balance. Thermal resistance model [27]. 


3.3. Conduction heat transfer through absorber tube 

Fourier’s law of conduction for a hollow cylinder is shown with 
the following equation [34] 

Q23cond = 27T/<23(T 2 — T 3 )/ ln(D 3 /D 2 ) (23) 


t?34 conv — 


2.425/c 34 (r 3 - r 4 )(PrRa D3 /(0.861 + Pr 34 )) 1/4 
(1 + (D 3 /D 4 ) 3/5 ) 5/4 


R0d3 


gjfjs - T 4 )D) 

otv 


The conduction coefficient is a function of the tube material. 
The current case comprises of three types of stainless steel, 304L, 
316L and 321H and a type of Copper B42. Conduction coefficient 
for 304L and 316L is [35] 

fc 23 = (0.013)T 23 + 15.2 (24) 



1 



3.6. Radiation heat transfer 


(26) 

(27) 

(28) 


for 321H 

k 23 = (0.0153)T 23 + 14.775 
and for B42 is 400 (W/mK). 

3.4. Heat transfer between absorber and glass envelope 



The radiation heat transfer between the absorber and the glass 
envelope is represented in the following equation 27,34]: 


Q34rad 


onD-i (t 4 3 - T 4 ) 

(l/e 3 + (l - s 4 )D 3 /(£ 4 D 4 )) 



The heat transfer between the absorber and the glass envelope 
is made up of radiation and convection. Convection heat transfer is 
a function of pressure [33]. When pressure is less than 1 mmHg, 
heat transfer is through molecular conduction unless, it is through 
free convection. 

3.4A. Convection heat transfer (vacuum) 

In vacuum state between the absorber and the glass envelope 
(pressure low than 1 mmHg), convection heat transfer is negligible 

[33], 

Q34 conv = 0 


3.7. Conduction heat transfer through the glass envelope 

For obtaining the conduction heat transfer through the glass 
envelope, the same equation (Eq. (23)) of the absorber are applied. 
It should be noticed that the temperature distribution is linear and 
the glass thermal coefficient is 1.04 [36]. 

3.8. Heat transfer between the glass envelope and environment 

Heat transfer between the glass envelope and environment in¬ 
cludes convection and radiation. The convection can be forced or 
natural depends on the wind condition. 


3.5. Convection heat transfer (with pressure) 

When pressure is over 1 mmHg, the convection heat transfer is 
gained with the following equations [27,32 . 


3.8A. Convection heat transfer 

Eq. (30) shows the convection heat transfer between the glass 
envelope and atmosphere [27] 

Q 56 conv = hsQDs7l(Ts — Ts) (30) 
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hse = Nu D 5 yl (31) 

E> 5 

Here, the Nusselt number is a function of wind condition. 


Case 1; without wind: 

Natural convection is observed when wind is not blowing. In 
this case, when 10 5 < Ra D5 < 10 12 , the following equations can be 
employed. 




0.38 7Ral 5 6 

[1 + (0.559/Pr 56 ) 9/16 ] 8/27 




g/HTs - r 6 )Ds 

$56^56 



solved and the mean temperature of the absorber tube is found. 
To find the glass envelope temperature, we used an iterative 
technique. First, we assumed a value for U L (initial guess for 
U L is 1 in our code) and then assumed a value for glass enve¬ 
lope temperature. An initial guess near the environmental tem¬ 
perature must be made for the glass envelope temperature and 
it was 30 °C as an initial guess. If the thermal difference be¬ 
tween the external surface of the absorber tube and the glass 
is not equal to the thermal difference between the glass surface 
and the environment, the procedure needs to go onto reach an 
equal thermal loss. Afterward, if the initial U L is different with 
U L through following equation, it has to be started from the 
first point. It is clear, for any assumed value of the U Ll a value 
for the mean temperature of absorber pipe is obtained by 
considering Eq. (18). 


P 





^56 

$56 



Case 2; with wind 

Forced convection is noticed when wind is blowing. In this case, 
when 0.7 <Pr 6 < 500 and 1 < Re D5 < 10 6 , the following equations 
are available (see Fable 1). 



CNu™Pr n 6 




Pr ^ 10 : n = 0.37 
Pr > 10 : n = 0.36 

All the fluid properties are taken into account at the atmo¬ 
spheric temperature, except Pr s which is at the temperature of 
the external surface of the glass envelope. 


3.8.2. Radiation heat transfer 

It is assumed that the cover is a small convex gray object in a 
large black body cavity, the sky. In this case, net radiation transfer 
between the glass envelope and sky is given by [27,37]: 

<j57rad = <X7lD 5 £ 5 (it - Tj) (37) 

To simplify the model, the effective sky temperature is approx¬ 
imated as 8 °C below ambient temperature. 


4. Result and discussion 

4.1. Thermal loss calculation 


TT (934 conv "T 934 rad) /oo\ 

Ul = 7rD 3 (r 3 -r 6 ) (38) 

A parabolic collector to generate power is modeled using MAT- 
LAB software with respect to optical and thermal analysis. A large 
number of characteristics including atmospheric pressure, envi¬ 
ronmental temperature and wind velocity can be employed in 
the presented model. With regard to required power and other 
restrictions, the absorber tube and the concentrator dimensions 
are accomplished. 


4.2. Concentrator 

In most of the collectors, the edge angle in concentrators is 
between 70 and 110 °C. The angle can be optimized with respect 
to manufacturing requirements and mechanical resistance [38]. 
In the present study, the edge angle is considered 90 °C which 
implies the minimized reflected beams and negligible tracking 
errors. The incident angle is normally assumed 0.95 for the men¬ 
tioned edge angle. 

Concentrator radiation coefficient is a function of the sur¬ 
face coating. Silvered glass and metalized acrylic are employed 
as the coating which have emissivity coefficient over 0.9. Met¬ 
alized acrylic is cheaper but has lower corrosion resistance 

[38]. 


4.3. Absorber tube 

The absorber tube is often made up of steel and covered by glass 
coating. Carbon steel absorber is well-known but 321H, 316L and 
304L stainless steels are used for high temperatures. Recently, 
321H has been used much since it decreases the bending problems 
and is resistant to hydrogen corrosion. 


The prominent point in thermal loss is that thermal loss 
from external surface of the absorber tube to the glass envelope 
has to be equal to the thermal loss from the glass envelope to 
the environment. In this regard, with initial fluid temperature 
and initial guess for U L , Eqs. (11), (12), (12a), (13)-(37) are 


Table 1 

Constants for Eq. (36). 


Rep 

C 

m 

1-40 

0.75 

0.4 

40-1000 

0.51 

0.5 

1000-200,000 

0.26 

0.6 

200,000-1,000,000 

0.076 

0.7 


4.4. Coating on the absorber 

The absorber is covered by a coating to maximize the absorp¬ 
tion coefficient and minimize the emissivity coefficient. 

Black chrome, black nickel and cermet are the popular coatings. 
Cermets have high temperature resistance. Black nickel has higher 
absorption coefficient and is cheaper but must be used less than 
300 °C. 

Four kinds of coatings are compared in Fig. 4 and their absorp¬ 
tion and emissivity coefficients are indicated in fable 2. 

The emissivity coefficient of black chrome is assumed as 0.11 
when it is less than 0.11. The emissivity coefficient of Luz-Cermet 
is 0.05 in a case that it is less than 0.05. Figs. 4-8 are gained using 
MATLAB software. 
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Black_Chrome 

— — Luz_Cermet 

■ (Solel_UVAC_Cermet(a 

• (Solel_UVAC_Cermet(b 



Ti (°C) 


Fig. 4. Efficiency as a function of fluid temperature and different coatings. 


4.5. Annular space pressure (between the absorber tube and the glass 
coating) 

Fig. 5 represents the collector efficiency as a function of 
pressure in the annular space between the absorber tube and the 
glass coating. To reduce the thermal loss, vacuum state in the 
annular space is a solution and heat transfer occurs only by 
radiation. 

As it is clear, the highest efficiency is the case of vacuum state in 
the annular space because there is not any convection heat trans¬ 
fer. The efficiency decreases with an increase in pressure. Because 
with increasing pressure, the heat transfer coefficient and conduc¬ 
tion coefficient and therefore heat loss, will be increased. Also, 
more pressure has more effect on the efficiency (the slope is 
greater). 

4.6. Gas in the annular space 

The collector efficiency for air, hydrogen, argon and vacuum 
state in the annular space is depicted in Fig. 6. The vacuum state 
has the highest efficiency. Inert gases are in the next grade. How¬ 
ever, is not much different with air. 

Filling the annular space with inert gas improves stability and 
reduces hydrogen leakage from the working fluid [33 . As shown 
in Fig. 6 slopes, respectively, for the case of vacuum mode, argon 
that is an inert gas, air and hydrogen increased, means their effect 
on efficiency, respectively, will be increased. Hydrogen has notice¬ 
able effect on the efficiency. 

4.7. Glass coating 

The glass coating reduces the thermal loss from the absorber. 
Low interval between the glass coating and the absorber causes 
high conduction thermal loss and high interval causes high convec- 


atm 1 — — atm 0.5 ■ atm 0.1 • vacuum * atm 0.01 



Fig. 5. Efficiency as a function of fluid temperature and annular space pressure. 



Fig. 6. Efficiency as a function of fluid temperature and the type of the gas. 


(m/s) 0.05 — — (m/s) 1 — >— (m/s) 3 • (m/s) 10 



Fig. 7. Efficiency as a function of fluid temperature and wind velocity. 


tion thermal loss. The optimal interval is a function of Rayleigh 
number but it is often considered 8 mm. 


Table 2 

Coatings absorption and emissivity coefficients. 


Coating Absorption coefficient Emissivity coefficient 


Black_Chrome 

0.94 

£ 3 = 0.0005333 x I 3 - 0.0856 


Luz_Cermet 

0.92 

£ 3 = 0.000327 x T 3 - 0.065971 


Solel_UVAC_Cermet (a) 

0.96 

s 3 = 2.249 x 10~ 7 x T\ + 1.039 x 10 4 

x I 3 +5.599 x 10 2 

Solel_UVAC_Cermet (b) 

0.95 

s 3 = 1.565 x 10 7 x T\ + 1.376 x 10 4 

x I 3 +6.966 x 10 2 
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Table 4 

Collector’s dimensions for different flows. 


Flow rate (kg/s) 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

0.1 

Collector length 

3.1 

3.1 

2.9 

3 

3.6 

3.1 

3 

(m) 

Collector aperture 

2.5 

2.4 

2.5 

2.4 

2 

2.3 

2.3 

(m) 

Collector power 

3502 

3499 

3495 

3499 

3498 

3501 

3504 

(W) 

Efficiency 

0.564 

0.587 

0.602 

0.607 

0.607 

0.607 

0.634 


8Ta= - 18Ta= ■ 27Ta= • 45Ta= 



Ti{°C) 


Fig. 8. Efficiency as a function of fluid temperature and environment temperature. 


Table 3 

Design parameters and operating conditions used for the case study. 


Direct normal incident solar irradiation (. I b ) (W/m 2 ) 800 

Name of fluid Therminol VP1 

Ambient pressure (atm) 1 

Wind speed (m/s) 3 

Ambient temperature (°C) 27 

Reflective 0.85 

Incident 0.95 

Edge angle (°) 90 


4.8. Environmental conditions 

Fig. 7 illustrates the efficiency as a function of wind velocity. 

It is obvious that the efficiency decreases with an increase in 
wind velocity. With increasing wind velocity, the effect of forced 
convection will be increased. Thus, it causes increasing heat loss. 
It is clear that there is no wind, there is only natural convection. 
Fig. 8 depicts the efficiency as a function of the environment 
temperature. 

It is inferred from Fig. 8 that the efficiency of the collector 
increases with increasing of environmental temperature. This 
increasing causes reducing the temperature difference between 
air and glass cover. Also, affects the Nusselt number and conse¬ 
quently convection coefficient. So increasing environment temper¬ 
ature decreases convection and radiation heat transfer from the 
glass cover to the environment. 

4.9. The variables and design factors 

The absorbed solar flux is a function of time. The program 
receives the time data (date). The program and is run for Tehran 
at 13:00 on July 2nd. The normal radiation is assumed at 800 W/ 
m 2 . The environment temperature, pressure and wind velocity 
are considered at 27 °C, 1 atm and 3 m/s, respectively. 

The working fluid inlet temperature is assumed at 100 °C and 
Terminol vpl is used as the working fluid. 

Stainless steel grade 321H is applied for the absorber tube and 
the outer diameter of the absorber tube is considered as 4 cm. 
Black chrome is employed as coating on the absorber tube. The 
annular space between the absorber and the glass coating is 
assumed 8 cm and filled with air at 1 atm. 

The reflective, incident and edge angles for the concentrator 
are assumed at 0.85°, 0.95° and 90°, respectively (see Table 3). 

Table 4 shows collector’s dimensions for different flows. 

In the present project, a parabolic collector is modeled using 
MATLAB software. Wind velocity, air temperature and other effec¬ 


tive parameters are discussed. The model is run for the maximum 
power at 3500 W. 


5. Conclusion 

Most collectors are designed to be coupled to a Stirling engine, 
are dish collectors. Also, they have high power or very low power 
and only one type of Stirling engine can couple to them. In this 
work, a collector is designed to connect to a Stirling engine with 
specified power. We obtain results such as collector length and 
aperture and other parameters for specified flow rate and condi¬ 
tions. Also, we obtain results of changing some conditions and 
parameters. Our collector has advantages. It is a parabolic-type 
and is designed for Stirling engine, in small scale and for domestic 
consumption with 3500 W of power. We can couple different type 
of Stirling engine to this collector. This system can use in other 
fields like air conditioning and solar water heater, in the following 
of generate electricity from Stirling engine. Most important thing, 
our parabolic collector has good and reasonable efficiency about 
sixty percent. This collector gives us many options to choose com¬ 
ponents and their material depending on financial issues and other 
things. 
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